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ABSTRACT 


NOAA-7 AVHRR data from April 1982 was used to perform global-scale analysis 
of acrosol particle characteristics. Mahonv’s (1991) water vapor correction was incor- 
porated into the AVHRR multichannel satellite data processing technique used bv Frost 
(1988). Channel 4’Channel 5 brightness temperature difference was used as an estimate 
of water vapor content in an air column. Greatest measured watcr vapor content was 
at the equator, decreasing toward the poles. Applying the correction reduced the aerosol 
particle size index, resulting 1n an increase 1n aerosol] optical depths. In addition, the high 
particle size index over the low latitudes in the southern hemisphere noted by Frost, was 
reduced significantly after applving the correction. A comparison was made between 
Frost's values for the variable two-term Henyey-Greenstcin phase function and those 
derived with the water vapor correction. The basic shape of the curve is similar to that 


determined by Frost. 
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I. INTRODUCTION 


A. PURPOSE 

The availability of satellite-measured radiances provides an effective and relatively 
inexpensive means for determining aerosol particle distribution on a global scale. Satel- 
lites provide extensive coverage of the earth’s surface with sufficient temporal periodicity 
to allow for the study of large-scale global effects while at the same time providing 
small-scale local measurements desirable for military applications. Data from 
meteorological satellites require careful processing in order to distinguish the much lower 
radiances associated with aerosol optical depth variations (Durkee 1984) from the 
greater radiances associated with cloud albedos, satellites were designed to measure. 

Improved aerosol optical depth retrieval algorithms have facilitated ongoing efforts 
to establish a global aerosol particle data base. Frost (1988) developed an algorithm to 
derive aerosol particle characteristics from the Advanced Verv Thigh Resolution 
Radiometer (AVHRR) on NOAA polar orbiting satellites. Aerosol particle size index 
can be derived by a comparison of the sensitivities of two adjacent spectral windows and 
so detect aerosol scattering. Frost examined the performance of a variable phase func- 
t10n as an improvement to the constant phase function used previously. Mahony (1991) 
examined the effects of water vapor on aerosol optical depth retrieval. An estimate of 
water vapor content in a column of air can be derived by comparing the brightness 
temperature of two AVIIRR windows in the thermal infrared region. He concludes that 
this effect can be compensated for, eluninating an error of 15% 1n the phase function. 

The purpose of this study is to incorporate Mahonv’s improvement into the algo- 
rithm used by Frost, which did not account for water vapor appreciably affecting re- 
trieved radiances. The April 1982 data examined by Frost will be re-evaluated to 
compensate for water vapor. The variable two-term Henvey-Greenstein phase function 
recommended by Frost will be examined in light of the new, water vapor corrected par- 


ticle size phase function and compared to the values he generated. 


B. MOTIVATION 

Aerosols consist of a variety of solid or liquid material dispersed in the atmosphere; 
examples include volcanic ash, smoke, and pollen. These particles vary in their distrib- 
ution and concentration in the atinosphere, which can impact local and global radiation 


budgets. The extreme concentrations of smoke particles resulting from the oil fires in 


Kuwait created a laver that reduced daylight to night in the region. Even modest con- 
centrations have demonstrated an impact in models of the earth’s radiation budget 
(Weiss et al. 1974; and Coakley 1976). A Soviet study of longwave radiation in the 
tropical Atlantic ocean (Zaitseva 1976) linked short-term temporal changes in temper- 
ature to changes in aerosol particles. 

Charlson et al. (1987) proposed that the major source of cloud condensation nuclei 
(CCN) over remote unpolluted oceans appears to be a waste product of phytoplankton, 
aqueous dimethylsulphide gas (DMS). A complex chain of chemical and biological 
processes result in DMS passing through the sea surface. Once 1n the atmosphere, it 1s 
oxidized to form sulfate aerosol particles which act as CCN. Varitions in phytoplankton 
populations with the associated changes in CCN particles could alter the earth’s heat 
budget. 

Of direct importance to the military is the impact of aerosols on electro-optical de- 
vices. The .4 to 1.0 wt waveband 1s considered the optimum window for lasers as it 
avoids ozone absorption at shorter wavelengths and water vapor absorption at longer 
wavelengths (Bloembergen et al. 1987). Aerosol scattering is a significant factor in this 


Same waveband. 


Il. THEORY 


A. DIFFUSE RADIANCE 
Knowing that optical depths are small over the ocean and that single scattering 1s 
dominant (Pfeil 1986), and assuming an optically thin atmosphere, Frost (1988) used the 


following simplified version of the Radiation Transfer Equation (RTE): 





ORLA 
L(é6, pu, d) ~ an p(O)d, pas 


where: 


Ie 


diffuse radiance, 


& 
I 


azimuth angle, 

Wy, = single scatter albedo, 

F, = solar irradiance, 

uw = cos0@,(0 = satellite zenith angle), 
p = scattering phase function, 

© = single scatter angle, 


6 = optical depth. 


Griggs (1975) and Durkee et al (1986) demonstrated this linear dependence of radiance 
on optical depth. 

There are three primary sources of satellite-detected radiance; aerosol scatter (L,), 
Rayleigh scatter (L,), and sea surface reflectance (L,). This relationship can be expressed 


as 
(Fee Tee 2.2 


Ravleigh scatter occurs when the size of the aerosol particle 1s much smaller than the 
wavelength of the incident radiation and, though strongly wavelength dependent. it does 
not Vary spatially. Sea surface reflectance is very small for red wavelengths, the albedo 


is 0.5%, and zero for radiation with wavelengths longer than 0.7m, except Where the 


sun-earth-satellite geometry results in sunghint (Ramsey 1968). Therefore, by avoiding 
data contaminated by sunghnt, total radiance can be approximated by the sum of the 
radiances due to Rayleigh and aerosol scatter (L = L, + Lz). By subtracting the compo- 
nent due to Rayleigh scatter, aerosol particle radiance (L,) will remain: 


a Wolo 
L,=L(e,u,@)— Lp~ ae P(O)d 4. 2 


where 0, 1s the aerosol particle optical depth. 

With Rayleigh scatter forming a baseline radiance, aerosol scatter is the dominant 
source of upwelled radiance variations ( AL ~ AL, ). Eqn 2.3 indicates satellite measured 
radiance is directly related to the changes in the aerosol particle characteristics of wa, 
p(O) and 6,. It 1s apparent that a knowledge of p(©) variations 1s required to determine 


a more,accurate o, [rom L,. 


B. OPTICAL DEPTH 
Optical depth is a measure of the scattering and absorption material present in the 
atmosphere. It is determined by vertically integrating the extinction coefficient (¢,,,), 


from the earth’s surface to the top of the atmosphere (11), expressed by 


H 
é -| Geena. 2.4 
0 


Extinction is due to both absorption (a,,,) and scattering (o,,,,) of radiant energy; how- 


scot 


ever, in the marine environment, absorption can be neglected, leaving extinction due to 


sCattefine alone la. = g..) suem that 


H 
0) -| Gren, Ga: 2 
0 


The amount of scattering is determined by the size of the particles, their scattering 


efficiency, and the size distribution of particles in the column given by 


oO 
—— | tr Ova (My 2,r)n(r)dr, 2.6 
0 


where, 


mr? = particle cross sectional area, 
O(n. 4,7) = scattering efficréency, 
m = complex index of refraction, 
2 = wavelength, 
r = particle radius, 


nr) = number distribution of particles. 


Aerosol particles are most efficient in scattering radiant energy at wavelengths near 
their radius, becoming Jess efficient as the difference between the two increases. Any 
changes in the size distribution will produce a corresponding change in 6, and im the 
observed radiance L,. The variations in radiance between two wavelengths, specifically 


red and near-infrared, can be quantified by the ratio (Durkee et al. 1991) 


= Coe = [orop(©)0 4 Tred 
Cyan [ogp(Q)0 4 Inir ) 


No 
~) 


S12 


This ratio (S,,) is called the particle size index and is sensitive to variations in the slope 
of the particle size distribution, increasing as the slope increases. Small particles have the 
greatest influence if the slope is steep, while large particles are most 1mportant for more 
Cradualsiopes. Generally, (4), , due to sinall particle scattering, 1s greater than {i5,),, 


resulting in an S,, ratio larger than | (Mahony 1991). 


C. WATER VAPOR ABSORPTION 

The dominant influence on retrieved radiances in the visible and near-infrared 
wavelengths 1s aerosol scattering. Though water vapor effects are not as significant as 
aerosol scatter, they can result in an over-estimation of S,, which could be interpreted 
as a change in aerosol characteristics or on the other hand, could mask a significant 
change. As the desired result is the estimation of aerosol particle characteristics from 
satellite observations, inclusion of water vapor effects will improve the accuracy of re- 
sults. 

In Figure 1, the effect of water vapor can be seen in the transmittance (tr) for 
AVHRR channels 1 and 2. Note that channel | t 1s greater than channel 2, indicating 
less absorption is taking place. Therefore, L, would tend to decrease more due to ab- 


sorption by water vapor than L, leading to the over estimation of S,,. 
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Fieure 1. Atmospheric transmittance due to water vapor: | Transmittance 1s 
shown with NOAA & AVIIRR channels 1 and 2 response functions 
(Mahony 1991). 


D. BRIGHTNESS TEMPERATURE DIFFERENCE 

An estimate of the water vapor content in an air colunin can be made using the 
AVIIRR 10 to 13 xm window region. In channels 4 and 5, at wavelengths of 11 and 12: 
um respectively, absorption by water vapor is dominant while aerosol scattering is rela- 
lively insignificant. The benefits of using these channels is the colocation and temporal 
linkage with channels 1 and 2, providing a water vapor snapshot along the same path 
as the aerosol optical depth (6,) measurements. Although not the ideal sensor windows 
for water vapor content, their use is justified within the accuracy of the approximations 


made (Mahony 1991). 


The brightness temperature (7) is determined from the radiance received, due to 
direct transfer, from the ocean surface to the satellite sensor. Therefore, the 
transmittance (rt) is directly proportional to the retneved 7. Only absorption reduces the 
radiance reaching the sensor; scattering 1s relatively insignificant. 

Dalu et al. (1981) proposed a method for determining sea surface temperature using 
briglitness temperature difference corrected for air column water vapor content. After 
manipulation of Dalu’s formula, the vertically integrated water vapor density (1) 1s 


provided by 


ta 
CO 


w= A(T, — Ts) cos @. 


where A is a constant. Using A = 19600g° A“ #177, a scanning angle 0 = 0, and a variety 
of temperature and humidity profiles, Dalu (1986) apphed this relationship to a radiative 
transfer model. The resulting line graph (Figure 2) has a correlation coefficient of 
R= 0.99, with an error of + L.S5kg¢ m-? in density. 

Mahony determined that the water vapor influence on the particle size index was a 
function of w and S,,(dry), the latter being the magnitude of S,, when w=0. The 


equation for the curve that best fit his data 1s 


Sy) 
Pr O82 


ta 
\O 


Sa = 


Where S,, and w are determined from measured radiances. 


FE. VARIABLE PHASE FUNCTION 

Based on Frost's work, a variable single scatter phase function was used, having 
been determined to be more accurate 1n its optical depth estimation. Ihgher ratios, which 
are due to smaller particles, have miore backscatter than the lower ratios, due to larger 
particles. This provides the basic shape necessary for the phase function needed to esti- 
mate optical depth. 

The two-term Henvev-Greenstein scattering phase function (Lenoble 1985) was 


used. 


Water vapour (kg/m?) 





0 2 3 
T, - T, (deg K) 


Figure 2. Total water vapor content as a function of brightness temperature differ- 
cnee (ial 1936) 


HW'(1 — 9?) (1 — IF) = g5) 


A) aro 
Ce Se, cos 0) (1 + 23 4g cos 0) 


Where IV’ is a weighting factor, g, and g, are asymmetry curve factors, and @ ts the scatter 


angle. All three of these factors are functions of S,, written as 


V2 (025 se) mall 
g, = 1.04 — (0.08333 x S,,), De 
gy = 1.2 — (0.58333 x S))). 2.13 


From the equations 2.11 through 2.13, the phase function’s dependence on S$), 1s appar- 
ent. ‘Phe marine-type phase function describes large particle distribution, Jow S,, ratios; 


While the rural-type phase function follows small particle distribution, high S,, ratios. 


As the ratios change, the factors change, resulting in a varving phase fuction from one 
similar to a marine phase function to one similar to a rural phase function. Figure 3 


shows some of the variable phase functions p(Q) calculated by Frost (1988) from the 


indicated S,, values. 
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Figure 3. Aerosol particle size index and associated phase functions: (lrost 1988). 
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Mi PROCEDURES 


A. DATA AND PROCESSING 

The data examined was collected by the NOAA 7 polar orbiting satellite whose av- 
erage altitude was 833 km. This sun-synchronous platform crossed the equator at ap- 
proximately 0230 and 1430 local standard time (north- and south-bound, respectively). 
Of the five spectral bands or channels on the AVHRR instrument; channels 1, 2, 4, and 
5 were used in measuring upwelling radiance. Table 1 provides the bandwidths for all the 
channels (Jensen 1986) plus a weighted average of the solar irradiance (/,) over the 


bandwidth used for channels | and 2 calculations. 


Table 1. NOAA 7 AVHRR CHANNEL BANDWIDTHS. 


Bandwidth Solar Irradiance 
ic ee ear ye 1543536 ml on 
0.725-1.10 por 103.0908 vill cri? 












3.55-3.93 211 


1O.5-11.3 cae 


Ples=12.35 ju 








The archives of the National Center for Atmospheric Research (NCAR) in Boulder, 
Colorado, provided AVHRR Global Arca Coverage (GAC) data for two periods in April 
1982; the 5th - 10th and the 20th - 25th. April 1982 data was analvzed for comparison 
with the work done earlier by Frost (1988). 

All the data was processed on the Interactive Environmental Digital Analysis Lab- 
oratory (IDEA LAB) located at the Naval Postgraduate School, Monterey, CA. The 


IDEA LAB produced enhanced images and graphs from the processed data. 


B. ANALYSIS 

The data analysis scheme was based on that used by Pfiel (1986) and Frost (1988). 
The water vapor algorithm used by Dalu (1986) was incorporated to account for ab- 
sorption by water vapor, as recommended by Mahony (1991). Also accounted for in the 
algorithm were changes to the solar irradiance due to variations in the sun-earth dis- 
tance. During preprocessing, pixels located at latitudes greater than 70° or those over 
land were discarded. 

To remove sunglint, an empirical method using viewing geometry was applied. This 
method, a weighted sum of relative zenith and azimuth angles for the sun and satellite, 
identified approximately the same contaminated areas as the Avian sunglint routine, 
Which 1s based on the Cox and Munk (1954) method. 

Ozone absorption effects occur in the red and near infrared wavelengths with the 
greatest absorption at a wavelength of 0.694 (Fleagle and Businger, 1980), which 1s in 
the channel I window. This effect was inadvertently compensated for twice by Frost. The 
error has been corrected for this study. 

Each pixel was then examined for clouds as only clear pixels were to be processed. 


Cloud removal was performed by the same method used by Frost: 


1. Test for bad pixels -- If channel 2 or channel 4 radiance values are less than Zero, 
after Ravleigh radiance was removed, then the pixel is classified as “bad” and not 
processed. 


2. Test for high clouds -- A pixel was considered to contain high clouds if its channel 
4 brightness temperature was less than 273°K. 


3. Test for low clouds -- A pixel was considered to contain low clouds if its channel 
2 albedo was greater than 40%. 


4. Preliminary test for clear pixels -- A channel 1/channel 2 albedo was calculated and 
if that ratio was greater than 1.5, the pixel was tentatively classified as clear. 


5. Final test for clear pixel -- The channel 2 albedo of the pixel tentatively identified 
as clear was compared with the four pixels above, below, to the right and to the 
left. All the pixels must fall with in five counts of each other for the pixel to be 
classified as clear, otherwise it was counted as partly cloudy. A count is equivalent 
to 0.1069% albedo (Kidwell 1986). 
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The clear pixels were then processed for a variety of parameters. The list is basically 
the same as Frost’s, the most significant difference is the inclusion of two parameters 


with the water vapor correction, the terms include the word “dry” in their name. 
Channel 1/Channel 2 ratio, 
Dry Channel 1,Channel 2 ratio, 
Aerosol optical depth, 
Dry Aerosol optical depth, 
Average phase function of Channels | and 2, 
Channel 4;Channel 5 brightness temperature differences, 
Number of clear pixels, 
Number of partially cloudy pixels, 


lr Dabameter. 


The P parameter is used to examine the phase function calulated bv Frost. 

The phase function p(@, S,,) is an important factor in the determination of aerosol 
optical depth. Any changes made in the phase function algorithm, such as the addition 
of a water vapor correction to S,,, will produce corresponding changes to the calculated 
optical depths. A comparison of this modified phase function to that determined by 
Frost should provide confirmation that the changes made produce physically sound re- 
sults and add to the knowledge of p(@, S,,). 

A new parameter was defined using the simplified RTE (2.1) that 1s calculated from 


radiance data and can be used to derive p(@) and 6,: 


Sul. 4 


P= ie 


= pl, Si2)o4. 3.1 


Using the fact that integrating the phase function over a sphere is unity, the following 


formirecan be Used to determine pe, S,.): 


[ov $,,)dQ = | pce, $15)5 dQ = | arac = 1, 3.2 


where A is equal to the inverse of the optical depth. It is this equation from which the 


phase function will be calculated for comparison to Frost’s curves. 
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All of the results were stored and averaged in 1° by 1° boxes; each data point is one 
box. The number of pixels varied from box to box, depending on the number of partially 


cloudy to cloudy pixels found in each square. 


[Ve See Swi 


A. INTRODUCTION 

Ongoing efforts to determine acrosol particle characteristics have resulted in a 
steady improvement in the algorithm. Water vapor effects were shown to be an impor- 
tant factor in aerosol optical depth retrieval but have not been previously applied to 
satellite data. In this study the latitudinal variation in water vapor content is initially 
examined to give an overview of where the effects could be the most significant. Particle 
distribution characteristics are examined next, both with and without the water vapor 
correction. The data having the correction will be refered to as “dry”. Both the corrected 
and uncorrected particle size index and acroso} optical depth will be examined and 
compared over the Pacific Occan followed by the latitude profiles of the averaged valucs. 
Finally, the shape of the P parameter curves are compared to Frost's phase function 


CUITVGS. 


B. COLUNIN WATER VAPOR 

Figures 4 and 5 show the average Water vapor content in an air column to provide 
a gencral overview of water vapor cffects in the atmosphere for the periods of April 5 - 
10 and 20 - 25, respectively. The general structure is that of a maximum in the vicinity 
of the equator with values then decreasing toward the poles.The drop is more rapid in 
the northern hemisphere than the southern hemisphere, the latter is transitioning from 
its summer season when temperatures are higher and more water vapor is present. The 
greatcr year-round air and sea surface temperatures at the equator allow more water 
Vapor to exist in a column of air. From this it can be seen that the greatest potential for 
error duc to Water vapor 1s in the equatorial region. An examination of the particle size 


index before and after applying the water vapor correction will be made next. 


LATITUDE PROFILE 
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Figure 4. Brightness temperature difference for April 5 - 10. 
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Figure 5. 
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C.. PARFICLE SIZEAINDEX 

Figures 6 and 7 depict the uncorrected S,, values across the Pacific Ocean for April 
5 - 10 and 20 - 25, respectively. A broad region of high values extends across the low 
latitudes in the southern hemisphere in both figures, indicating the presense of a greater 
number of smaller particles. Low values indicating larger particles can be seen in Figure 
6 across the northern hemisphere. In Figure 7 there are two noticable maxima; one in 
the vicinity of the New Hebrides islands, the second west of South America. The north 
Atlantic appears to be a region of small particles except along the east coast of the 
United States. 

In Figures 6 and 7, bands cross each figure from north-northwest to south- 
southeast. These are a result of the moving of the sensor’s aperature to the left and right 
of the satellite’s path to view the side regions. Aerosol scatter 1s anisotropic, favoring 
forward and backward scatter. A satellite is looking at backscatter when it looks directly 
at the earth. But as the sensor looks to the side, it is no longer viewing backscatter, but 
the smaller scatter associated with smaller scattering angles. This results in a decrease 
in the phase function and a corresponding increase in the optical depth. Future work 
will need to comect for thus effect. 

Figures § and 9 show a comparison of the particle size index profiles for S,, and 
S,,(dry') for early and late April, respectively. The S,, profiles in both figures show a 
broad region of high values, approximately 2.35, from the equator to 40°S; the values 
then decrease toward the poles more rapidlv in the southern henusphere than in the 
northern. The greatest drop can be seen in the southern hemisphere in Figure 8 with a 
minimum value of 2.0. The minimum in Figure 9 1s 2.10. Frost’s values were 2.0 at the 
equator decreasing to 1.65 moving southward toward the pole; a second smaller decrease 


occurred toward the north pole, dropping to 1.9. His underestimation was due to a 
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Figure 6 Particle size index over the Pacific, April 5 - 10 
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Figure 9.) Comparison of particle size index profiles for April 20 - 25: the solid line 
is S,, and the dashed line ts $,,(ry). 
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In early April, Figure 10, distinct high value areas remain in the Gulf of Carpentaria 
and along the coast of China. Late April (Figure 11) shows two areas of higher values 
in the northern Atlantic Ocean and along the west coast of North America. The maxi- 
mum west of South America in Figure 7 is easier to distinguish in Figure 11 and has 
lower values. The maxima off the Americas appear to be associated with upwelling re- 
gions that are biologically very productive, which could be a source for the greater 
number of small aerosol particles. 

Mei ures o and Jit 1s easy 10 see the decrease tn size index values Between S,, and 
S,(dry). At the equator S,,(dry) 1s 0.5 less than S,, 1n early Apnil, 0.4 less later in the 
month. From the equator to the poles, the difference between S,, and S,,(dry;) decreases 
to approximately half the values at the equator in both figures. In both instances the 
overall shape of the S,,(dry) curve shows a minimum near the equator gradually in- 
creasing as you move toward the poles, which is the reverse of Frost’s results where 
values decreased toward the poles. Though a more uniform profile was expected, the 
peaks that exist 1n the corrected profiles of both figures correspond to the latitudes of 


the maxima depicted in Figures 10 and 11. 


FE. OPTICAL DEPTH 

The views of the Pacific in Figures 12 and 13, show the southern henusphere with 
low optical depths and the northern hemisphere with very high optical depths. The 
largest feature is the wide area of high values, 0.45 to 0.50, that occurs along the 
Mexican/Central American west coast extending westward into the Pacific. This feature 
isthe result ol the El Clicon volcame (17.3 495.27) \’ertpting on a nuniber of occa- 
sions during the period of 28 March to 6 April. The high values in the Gulf of Mexico 
@tee als tle result. ol “El Chicon./-Ihe mioveinent tothe east 1s possibly due to 


Stratospheric winds carrying particles ejected to those heights (Frost 1988). In early 
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Figure 10. Dry particle size index over the Pacific, April 5 - 10 


April (Figure 12) the plume extends only half way across the basin, later in the month 
It appears to approach Japan (Figure 13) 

A second feature can be seen over the East China Sea in Figure 13. Thus reg: 
high optical depth corresponds to a region of large particles in Figure 7, Gobi di 


blowing off mainland China produced this feature, as noted by Frost (1988) 
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Figure li. Dry particle size index over the Pacific, April 20 - 25. 
y | { 


Figures 14 and [5 show a comparison between the optical depth profiles for 0, and 
O,(dry), The 0, profiles both show a maximum value of 0.2 at approximately 20° 
which is due to the eruption of the EI Chicon volcano. A second peak located in 
high southern latitudes can also be seen. The reason for this maximum is un 


possibly the result of sea salt from the circumpolar winds around Antarcti 
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Fieure 12. Optical depth over the Pacific for April S - 10. 


KF. DRY OPTICAL DEPTH 

Figures 16 and 17 depict 6,(dry) across the Pacific Ocean. The high optical depth 
across the northern hemisphere due to El Chicon and Gobi dust have greater a: 
tent and higher values than those in Figures 12 and !3. These features cover regi hat 
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Figure 13. Optical depth over the Pacific for April 20 - 25 


and LS al ¢ LHe ry yptical de} th i } [1] Pip OAaFrYV) f 1g tne Carl ALLA j | \{ 
They show that the water v: por correction resul Pitan Hicre 
ues. Increases were generally on the 


slighth ereater changes occurring at preater aistances [ron 











LATITUDE PROFILE 


30° N 


10°S 


20-5 


30°S 


40°S 


SoS 


60°S 





ies 
0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 


Optical Depth 


Figure 14. Comparison of optical depth values for April 5 - 10: solid lines ure 0, 
Values, dashed are 6,(céry’) values. 


28 





LATITUDE PROF TLE 


70°N 





70°S 
0.00 0.05 0.10 0.15 0.20 0.25 0.30 90.35 0.40 


Optical Depth 


Figure 15.) Comparison of optical depth values for April 20 - 25: solid lines are 6, 
values, dashed are 0,(dry) values. 


changes occur in areas where the number of aerosol particles are sufficient to affect the 


optical depth. 


G. PHASE FUNCTION 

Figure 18 shows the P parameter as a function of scatter angle for three different 
particle size indices. A comparison to Frost’s plot of the Henyey-Greenstein phase 
function (figure 3) indicates that the basic shape of the P parameter curves are similar 
to that shown by Frost. A smaller particle size index has a larger value in the forward 
direction (between 0° and 60° ) than larger index values. The larger the particle size index 
the closer the curve approaches that of Rayleigh scatter, which 1s expected as smaller 
particles dominate. Where the scatter angle is between 90° to 140°, the P parameter ap- 
proaches a minimum value for each curve. As the angle increases to 180° the curves 


show an increase in P values. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


Frost showed that his technique for processing AVHRR data was effective as a tool 
in the analysis of aerosol particle characteristics. Though significant improvements were 
made with his study, sources of error have since been identified requiring futher im- 
provements to his algorithm. The determination that water vapor was a significant 
source of error suggested a re-evaluation of the April 1982 data in order to assess water 
vapor impact on aerosol optical depth retrieval. 

The results of this study have shown that inclusion of a water vapor correction 
produces a change of approximately 15% 1n the aerosol particle size index. Without this 
correction, large regions of high particle size index values, which correspond to a 
predominance of relatively smaller particles, were actually regions of lower values 
masked by water vapor effects. This made it difficult to see true areas of small particles 
as they were indistinguishable from the water vapor variations. 

Since the phase function algorithm 1s a function of S,,, the modification made for 
water vapor will affect the phase function. It has been shown that the shape of the 
modified phase function is similar to the one determined by Frost. The improved 
parameterization cannot be directly compared to Frost’s until the integration (3.2) 1s 
accomplished. 

The aerosol optical depth has also been shown to be modified by the water vapor 
correction producing an approximate change of 5%. The observed extent and magni- 
tude of areas of high dry optical depths was greater than the optical depths without the 
correction. 

Additional improvements to the algorithm used are still necessary. Multiple scatter 


has not been fully accounted for. The assumption used to allow for single scatter only 
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is valid for optical depths less than about 0.2. This value was exceeded in the regions of 
high optical depth associated with major aerosol particle producing events such as the 
E] Chicon eruptions. Errors mn the values due to multiple scatter exist and must be cor- 
Fecred, 

A second area for improvement would be in the determination of S,,. Current pro- 
cedures use a ratio of channel I’channel 2 radiances to calculate the particle size index. 
The preferred method is to use a ratio of channel I/channel 2 optical depths, which re- 
quires a better knowledge of the phase function for both channels. It 1s the optical depth 
which contains information on the aerosol particles and would provide the best estimate 
for S,,. In this study the phase function and optical depths are at the channel | wave- 
length. This improvement should eliminate the banding that occurs in the Pacific basin 
photographs which possibly hides features of importance. 

Finally, additional studies of global aerosol particle characteristics are needed. This 
study was limited to a single month. Studies exploring the changes that occur during a 
year and those consisting of a number of vears by month or season would increase our 


understanding of the role of aerosol particles in the global clrmiate. 
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